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ABSTRACT 

We use a new approach to obtain limits on the absorbing columns towards an initial sample 
of 10 long Gamma-Ray Bursts observed with BeppoSAX and selected on the basis of their good 
optical and nIR coverage, from simultaneous fits to nIR, optical and X-ray afterglow data, 
in count space and including the effects of metallicity. In no cases is a MW-like extinction 
preferred, when testing MW, LMC and SMC extinction laws. The 2175A bump would in 
principle be detectable in all these afterglows, but is not present in the data. An SMC-like 
gas-to-dust ratio or lower value can be ruled out for 4 of the hosts analysed here (assuming SMC 
metallicity and extinction law) whilst the remainder of the sample have too large an error to 
discriminate. We provide a more accurate estimate of the line-of-sight extinction and improve 
upon the uncertainties for the majority of the extinction measurements made in previous studies 
of this sample. We discuss this method to determine extinction values in comparison with the 
most commonly employed existing methods. 


Subject headings: galaxies: ISM — gamma rays: bursts 


1. Introduction 

The accurate localisation of Gamma-Ray Bursts 
(GRBs) through their optical and X-ray after- 
glows has enabled detailed studies of their environ- 
ments. Selection by the unobscured gamma-ray 
flash alone has allowed the discovery of a unique 
sample of galaxies which span an enormously wide 
range of redshifts from z ~ 0.009 (GRB 980425, 
e.g. Tinney et al. 1998; Galama et al. 1998a) to 
6.3 (GRB 050904, Kawai et al. 2005, 2006). The 
subset of long-duration (> 2 s) GRBs are almost 
certainly caused by the collapse of certain massive 
stars to black holes, confirmed by observations of 
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supernova components in the late-time afterglows 
of a number of long-GRBs (Woosley & Bloom 
2006; Kaneko et al. 2006) and by the observed lo- 
cation of GRBs in UV-bright regions within their 
host galaxies (Bloom et al. 2002; Fruchter et al. 
2006). GRBs are located in host galaxies which 
are generally small, faint, blue and highly star- 
forming (e.g. Chary et al. 2002; Fruchter et al. 
1999; Le Floc’h et al. 2003). Hence, detailed 
and extensive host galaxy observations provide a 
wealth of information on the gas and dust prop- 
erties of star-forming galaxies throughout cosmo- 
logical history. 

Accurately measuring the dust content of these 
galaxies is of great importance in, to name one ex- 
ample, the determination of their unobscured star 
formation rates where uncertainty in the correc- 
tion for dust can easily dominate the errors on 
the measured star formation rates for high red- 
shift galaxies (e.g. Pettini et al. 1998; Meurer 
et al. 1999). Absorption within our own Galaxy 
along a particular line of sight can be estimated 
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and removed, but absorption which is intrinsic to 
the GRB host galaxy as a function of wavelength 
is unknown, and is especially difficult to determine 
given its dependence on metallicity and the possi- 
ble existence of dusty intervening systems whose 
extinction curves cannot be disentagled from those 
of the host galaxy. Afterglow spectroscopy and/or 
photometry can be used to provide an estimate of 
the total extinction along the line-of-sight to the 
GRB. If the host galaxy itself is bright and ex- 
tended enough to be observed once the afterglow 
has faded, different lines-of-sight may be probed 
besides the off-centre UV-bright regions within 
which GRBs are generally situated. 

Extinction in the optical/UV regime due to 
dust grains is typically modelled using either 
Milky Way (MW or Galactic), Large Magellanic 
Cloud (LMC) or Small Magellanic Cloud (SMC) 
extinction curves (e.g. Pei 1992) because these 
curves can be measured and so are well known, or 
with the Calzetti extinction law derived empiri- 
cally from UV observations of starburst and blue 
compact galaxies (Calzetti et al. 1994). It appears 
that dust content in most GRB hosts produces an 
SMC-like extinction law (e.g. Galarna & Wijers 
2001; Vreeswijk et al. 2004; Stratta et al. 2004; 
Kann et al. 2006; Schady et al. in preparation), 
owing to an observed lack of the 2175A feature 
thought to be caused by carbonaceous dust grains 
(Draine & Lee 1984). This feature has, however, 
been clearly observed in GRB afterglow spectra 
where the line of sight between us and the GRB 
is intercepted by intervening systems: the best 
example to date, in which the extinction curve 
of an intervening system could for the first time 
be disentangled from that of the host galaxy, is 
GRB 060418 (Ellison et al. 2006). 

In general, low amounts of optical extinction 
are found towards GRBs, unexpected if GRBs 
are located in dusty star-forming regions, whilst 
the X-ray spectra reveal a different picture (first 
noted by Galarna & Wijers 2001). At X-ray wave- 
lengths absorption is caused by metals in both gas 
and solid phase, predominantly oxygen and car- 
bon (see e.g. Wilms et al. 2000), and we often 
measure high values for the absorbtion columns. 
However, these absorption edges are shifted out 
of the X-ray observing window for high redshifts 
(z > 2) beyond which only large columns can be 
measured and there is a degeneracy between red- 


shift and X-ray column density (e.g. Watson et 
al. 2002). The GRB host metallicity is observed 
to be low (compared with the Milky Way) in mea- 
surements via optical spectroscopy of a dozen or 
so afterglows. Host metallicities can reach values 
as low as 1/100 Solar (GRB 050730, Starling et 
al. 2005; Chen et al. 2005, but see also Prochaska 
2006 for potential caveats) - even lower than found 
for the SMC (see Figure 3 of Fynbo et al. 2006 for 
an overview). This only increases any measured 
X-ray column, which is expressed as an equivalent 
hydrogen column density, Ah- But here we note 
that metallicities are not generally obtainable for 
lower redshift GRBs (z < 2) due to the hydrogen 
Lyman-a line lying in the far UV outside typical 
observing windows. 

The apparent discrepancy between optical and 
X-ray extinction resulting in high gas-to-dust ra- 
tios in GRB host galaxies (often far higher than for 
the MW, LMC or SMC, e.g. GRB 020124, Hjorth 
et al. 2003, but see Schady et al. in prepara- 
tion) is not satisfactorily explained, though the 
suggestion that dust destruction can occur via the 
high energy radiation of the GRB (e.g. Waxman 
& Draine 2000) could possibly account for the dis- 
crepancy. It is thought that circumburst dust may 
be destroyed by sublimation of dust grains due to 
UV emission (Waxman & Draine 2000), sputter- 
ing (Draine & Salpeter 1979) or dust grain heat- 
ing and charging (Fruchter et al. 2001). Alter- 
native models for the extinction by dust grains, 
including skewing the dust grain size distribution 
towards larger grains, have been investigated, and 
in fact such a grain size distribution may result 
from exposure of the dust to the GRB radiation 
field since destruction of small grains is more ef- 
ficient than for larger grains (Perna et al. 2003). 
Attempts to model the process of dust destruction 
have been made by e.g. Perna & Lazzatti (2002). 
However, such models have not replaced SMC-like 
(low metallicity) extinction as the best descrip- 
tion of most GRB environments (e.g. Stratta et 
al. 2004; Kann et al. 2006). 

Traditionally the optical and X-ray spectra 
have been treated seperately in extinction stud- 
ies. Since the underlying spectrum is likely a 
synchrotron spectrum (power law (pi) or broken 
power law (bknpl), e.g. Galarna & Wijers 1998) 
extending through both wavelength regions, it 
would be most accurate to perform simultaneous 
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fits. More recently such fits have been made, ei- 
ther by fitting the X-ray spectrum individually 
and thereby transforming the model counts to 
flux to create a spectral energy distribution (SED) 
with the optical data in flux space (e.g. Stratta 
et al. 2005; Watson et al. 2006), or by using 
the Swift XRT X-ray and UVOT Ultraviolet (UV) 
and U,B,V band data together in a fit to the 
count spectrum (e.g. Blustin et al. 2005; Schady 
et al. in preparation). The Swift UVOT data 
can be loaded directly into the xspec spectral fit- 
ting package and treated in count space owing 
to its calibration, which is generally not true for 
ground-based data. We present here an alterna- 
tive method, which makes use of simultaneous fits 
in count space extending from near-infrared (nIR, 
in this case K band) to X-ray (10 keV) to ob- 
tain the most accurate possible measurements of 
both the underlying continuum spectrum and the 
extinction. 

In Section 4 we fit the broad-band SEDs (from 
nIR through X-ray) of a subsample of the Bep- 
poSAX sample of GRB afterglows to better mea- 
sure the extinction properties of their host galax- 
ies - a sample chosen for its availability of suitable 
data well studied in the seperate optical and X-ray 
band passes. Section 2 outlines the data sample 
and reduction techniques. Section 3 describes the 
method we use to model the broad-band SEDs 
and Section 4 presents the results of our fitting 
through discussion of individual bursts and com- 
parison with previous studies. In Section 5 we 
discuss the implications of our findings for galaxy 
extinction curves, and summarise and compare the 
various methods now available to measure extinc- 
tion in the hosts. We conclude by summarising 
our method and findings in Section 6. An analysis 
of the blastwave parameters and density profiles 
for the circumburst media obtained from these fits 
will be presented in a forthcoming paper (Starling 
et al., Paper II). 

2. Observations 

This sample of 10 long GRBs observed with 
the BeppoSAX Narrow Field Instruments (approx. 
0.1-10 keV) is chosen for the good availability (3 
bands or more) of optical/nIR photometry (Ta- 
bles 1 and 2). The optical/nIR bands available for 
each source and their references are listed in Ta- 


ble 3. As these GRBs are all previously studied, 
overlapping with the samples studied by Galama 
& Wijers (2601) and Stratta et al. (2004), this 
constitutes a good sample on which to first adopt 
this method of simultaneous SED fitting. 

All X-ray observations are taken from the Bep- 
poSAX data archives, using the LECS and MECS 
instruments raw data within the energy ranges 
0.1-4 and 1.0-10 keV respectively. Data have been 
reduced using the SAXDAS routines. We com- 
bined data from the MECS2 and MECS3 instru- 
ments (except in the case of GRB 970228, where 
we use the MECS3 instrument only, see Stratta et 
al. 2004), including a gain equalisation. We then 
combined multiple observations for each source 
and instrument type among the narrow field in- 
struments, omitting the last observation if it was > 
3 days later than the previous one, before extract- 
ing spectra. Background X-ray spectra were taken 
from blank fields and count rates checked against 
the local background finding no adjustments nec- 
essary: the net count rates of the two types of 
field agree on average to within 0.0001 counts s~ 1 . 
The latest canned arf and rmf files were used with 
MECS data, but were created for LECS observa- 
tions at their off-axis source positions (listed in 
Table 3 of Stratta et al. 2004). We group all spec- 
tra such that a minimum of 20 counts are in each 
bin in order to use the y 2 statistic. However, in 
the case of GRB 000926 there were very few counts 
in the X-ray spectrum so we have required only 10 
counts per bin. There is a known offset between 
the normalisations of the LECS and MECS instru- 
ments. We fit for this offset in the X-ray spectra 
only, adding a constant-value free parameter to 
the model and adopting MECS as the reference 
for the LECS instrument. We fix the offset values 
in the SED of each GRB to these values (Table 1). 

All temporal decay slopes, both for X-ray and 
optical lightcurves, have been taken from the liter- 
ature and are listed together with Galactic extinc- 
tion corrections in Table 2. Optical and nIR pho- 
tometry was taken from the literature and from 
our own nIR observations of GRB 990510 (de- 
scribed in Curran et al. in preparation). 

3. Method 

Per source all data are fitted simultaneously, 
assuming wherever possible no prior model. This 
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Table 1: BeppoSAX data table for combined datasets. We have measured the MECS-LECS offset values 
from the combined X-ray spectra with respect to the MECS instrument. 


GRB 

obs start obs end 
(days since 

obs midpt(log) ton-source (s) 

trigger) MECS LECS 

MECS-LECS offset 

970228 

0.344 

0.693 

0.520 

1.43x 10 4 

5.611xl0 3 

1 4K+U-67 

970508 

0.434 

3.091 

1.679 

5.90xl0 4 

2.36xl0 4 

i qq+0.50 

l.UJ— o 43 

971214 

0.274 

2.528 

1.362 

l.OlxlO 5 

4.62 xlO 4 

i 57+0.67 
1 *° ' -0.56 

980329 

0.294 

2.026 

1.148 

6.38xl0 4 

2.49 xlO 4 

0-68+3Q 

980519 

0.406 

1.468 

0.930 

7.82xl0 4 

2.31 xlO 4 

1-00+;! 4 

980703 

0.827 

1.902 

1.333 

3.92 xlO 4 

1.66xl0 4 

0-99+it 

990123 

0.242 

2.573 

1.245 

8.20xl0 4 

2.80xl0 4 

0.77±0.9 

990510 

0.334 

1.850 

1.067 

6.79xl0 4 

3.17xl0 4 

0-86+12 

000926 

2.003 

2.466 

2.234 

1.96X10 4 

5.027xl0 3 

0.70+ 43 

010222 

0.376 

2.703 

1.511 

8.84 xlO 4 

5.11xl0 4 

1-43+11 


is achieved by fitting in count space (as is tradi- 
tional in the X-ray regime where one fits for the 
emission model, extinction and instrumental re- 
sponse simultaneously): the optical and nIR mag- 
nitudes are converted to flux and then to counts. 
For the magnitude to flux conversion we use the 
zero points and effective bandwidths of each opti- 
cal band (Johnson for U, B, V, R, I, J, K, 2MASS 
for H and K s and Bessel for V c , R c , I c ). In the 
small number of cases for which the specific band 
is not stated, we assume the appropriate Johnson 
filter. These fluxes are then converted to photons 
cm~ 2 s -1 per bin (bin width = effective band- 
width of the filter) within the ISIS spectral fitting 
program (Houck & Denicola 2000) which is equiv- 
alent to the X-ray units of counts cm -2 s -1 per 
bin since total number of counts is conserved. 

Since we fit in count space, we need not first 
assume a model for the X-ray spectrum to convert 
the counts to flux. Herein lies one advantage of 
using this method. The second advantage comes 
through the multiwavelength approach. The nIR, 
optical and X-ray spectra are related since we as- 
sume the broadband spectrum is caused by syn- 
chrotron emission, hence a simultaneous fit pro- 
vides greater accuracy and consistency between 
the parameters. Inclusion of nIR data and R band 
optical data together with the 2-10 keV X-ray 
data, regions over which extinction has the least 
effect, allows the underlying power law slope to be 
most accurately determined. 

The X-ray data typically have much longer ex- 


posure times than the individual optical/IR mea- 
surements, particularly given that where sensible 
we have combined X-ray datasets from different 
epochs to increase the signal to noise. We have 
chosen to omit the last X-ray observation of any 
GRB which occurs more than 3 days after the pre- 
vious observation where this contributes more to 
the background noise than the signal and skews 
the observation time midpoint. We note that 
structure in the X-ray lightcurves and spectral 
changes at early times have often been reported 
for the better and earlier sampled Swift GRB af- 
terglows (e.g. Nousek et al. 2006), which may 
affect the BeppoSAX afterglows with fairly early 
coverage, i.e. those beginning at 0.2-0. 3 days since 
burst (see Table 1). All optical data are scaled to a 
common time, which corresponds to the midpoint 
of the X-ray observations, calculated in log-space. 
This is done by extrapolation according to the de- 
cay rates and lightcurve break times determined 
in a thorough analysis of pre-Swift burst optical 
lightcurves by Zeh et al. (2006, Table 2). 

We fit for the X-ray and optical extinction at 
the redshift of the GRB. For GRBs 980329 and 
980519 the redshift is not known, and therefore ab- 
sorption at the source cannot be measured. How- 
ever, a photometric redshift has been made for 
GRB 980329 of z ~ 3.6 (Jaunsen et al. 2003) 
which we adopt here. For GRB 980519 we adopt 
the mean redshift for this sample, z ~ 1.58, to 
make an estimate of the intrinsic absorption re- 
quired. 
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Flux is depleted in the bluemost optical bands 
for high redshift bursts due to the Lyman absorp- 
tion edges and the resulting transmission is calcu- 
lated for each band per burst (see Curran et al. 
in preparation for further details of this calcula- 
tion, which we note involves assuming a spectral 
slope for the optical flux over a. given band and 
uses parameters from Madau 1995 and Madau et 
al. 1996). However, this has an effect only on the 
U through R band magnitudes of GRBs 971214 
and 980329 and a minor effect on the U band 
magnitudes of GRBs 980519, 990123, 990510 and 
000926. 

The errors on the optical magnitudes are taken 
from the photometric errors in the literature, and 
are set to 0.1 for cases where the literature reports 
a smaller error to account for systematic uncer- 
tainties. The conversion from magnitudes to flux 
has an associated error of up to 5 % (Fukugita et 
al. 1995). We note that by extrapolating mag- 
nitudes to different times we introduce a possible 
random error, since there is an uncertainty in the 
measured decay indices and break times, which we 
can allow for in an offset between the optical and 
X-ray data. This error is not included in the in- 
dividual data points because this would introduce 
an artificially large error in the optical slope which 
is not in fact present since the relative errors be- 
tween the bands do not change. 

3.1. Models 

Data are fitted within the spectral fitting pack- 
age ISIS (Houck & Denicola 2000) using both mod- 
els written for use within xspec (Arnaud et al. 
1996) and models written for ISIS. Models consist 
of either a single or a broken power law, to allow 
for a possible cooling break in between the optical 
and X-rays. Should the break in the power law be 
due to cooling, the difference in slope is A/3 = 0.5 
(e.g. Wijers & Galama 1999), which we fix in the 
broken power law model. 

On its way to us, the intrinsic power law is ab- 
sorbed by optical extinction at the host/burst red- 
shift (this could be local to the GRB itself or an- 
other location within the host galaxy). The ex- 
tinction curves used for intrinsic optical extinction 
in this study are Galactic-, SMC- and LMC-like 
(Figure 1) following the prescriptions of Cardelli 
et al. (1989) and Pei (1992) respectively. We 
do not use the Calzetti extinction curve (Calzetti 


et al. 1994) because it has a larger error associ- 
ated with it, being constructed from fewer mea- 
surements than 'chose for the nearby Magellanic 
Clouds. There is also absorption in the X-ray 
regime predominantly by metals, e.g. the oxy- 
gen edge, for which we use a photoelectric ab- 
sorption model. We refer the reader to Wilms et 
al. (2000) for a detailed description of the X-ray 
absorption models. The X-ray absorption model 
can be computed for various metallicities by sim- 
ple scaling of the Solar abundances by a constant 
factor. We adopt firstly Solar metallicity and sec- 
ondly the metallicity assumed in the optical ex- 
tinction model: using SMC-like absorption one 
would adopt Z— 1/8 Z© and for LMC-like absorp- 
tion Z= 1/3 Z© (Pei 1992), for self-consistency. 

The flux is then corrected for Galactic absorp- 
tion (Table 1). In the X-ray regime these val- 
ues are fixed at the Nh values given in Dickey 
& Lockman (1990), which are averages over 1 
degree at the positions as given in the Simbad 
catalogues or from the BeppoSAX Narrow-field 
instruments. For completeness and comparison, 
we also list in Table 1 extinction measurements 
from the newer and slightly higher resolution Lei- 
den/Argentine/Bonn Galactic H I Survey (resolu- 
tion of ~0.6°, Kalberla et al. 2005). These new 
values are not significantly different than those 
of Dickey & Lockman though appear to be gen- 
erally lower, and to date all previous studies use 
the values from Dickey & Lockman which we will 
also use here. For the optical extinction we use 
E(B — H)Gai values given by Schlegel et al. (1998) 
from their full-sky 100 pm map together with the 
Galactic extinction curve of Cardelli et al. (1989) 
with R v = A V /E(B - V) = 3.1. The Schlegel et 
al. (1998) maps have a best resolution of 6.1’: for 
each source we use the best resolution available, 
but in some cases we must use an average over 1 
square degree centred on the source coordinates. 
The fit statistic calculated is x 2 , using a Levenberg- 
Marquardt fit minimisation method. Errors in the 
LECS-MECS offsets (except for GRB 970508), op- 
tical decay slopes and redshifts are not propagated 
through the fitting routines. These values are sim- 
ply fixed from Tables 1 and 2 respectively. We also 
do not include uncertainties on the zero points on 
our photometric data points (see previous Sec- 
tion). 
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Table 4: Results (main parameters) of fits to the spectral energy distributions. For each GRB we fitted the SED with the power law and 
the broken power law models with all 3 extinction models (MW, LMC and SMC). X-ray column densities given in brackets are for the 
appropriate metallicities (LMC or SMC), otherwise Solar metallicity is assumed. The LECS-MECS offsets were fixed at the central values 
listed in Table 2. For the broken power law models Id = r 2 - 0.5. Where the break energies are unconstrained (indicated by the letter U) 
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Table 4: - continued 
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Fig. 2. — Data (black crosses) and best fitting models (red lines, see Table 4) for each of the GRBs in the 
sample. Data (nIR and optical photometry and BeppoSAX LECS and MECS X-ray spectra) and models 
are shown in count space. The bin size (effective bandwidth) of the optical data points can be seen in the 
model fits. The lower panels show the deviation from the model for individual data points, in units of their 
measurement error. 
















Table 3: Range of optical and nIR magnitudes used 
in this study. f 


GRB 

bands used 

refs 

970228 

VR C I C 

[1] 

970508 

UBVR C I C K S 

[2,3] 

971214 

VRIJK S 

[4, 5, 6, 7] 

980329 

RIJK 

[8] 

980519 

UBVRJc 

[9] 

980703 

RIJHK 

[10] 

990123 

UBVR1HK 

[11] 

990510 

BVRIJHKs 

[12,13] 

000926 

UBVRIJHK 

[14] 

010222 

UBVRIJ 

[15] 


1 Galama et al. 2000 2 Galama et al. 
1998b 3 Chary et al. 1998 4 Halpern et al. 
1998 5 Diercks et al. 1997 6 Tanvir et al. 
1997 7 Ramaprakash et al. 1998 8 Reichart 
et al. 1999 and references therein (their Ta- 
ble 1) 9 Jaunsen et al. 2001 10 Vreeswijk et 
al. 1999 11 Galama et al. 1999 12 Stanek et 
al. 1999 13 Curran et al. in prep. 14 Fynbo 
et al. 2001 15 Masetti et al. 2001 
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Fig. 1. — The well known extinction curves for the 
Milky Way (MW), Large and Small Magellanic 
Clouds (LMC and SMC respectively, Pei 1992). 
The transmission of optical/UV intrinsic flux with 
energy is shown for an object at redshift z = 0 and 
with a large optical extinction of E(B — V) = 0.2. 


4. Results 

The SEDs and the results of fits t® the SEDs for 
all GRBs in the sample are listed in Table 4, and 
best fitting models are shown overlaid on the data 
in Figure 2. Figure 3 shows a comparison of the 
absorption measurements with Galactic, LMC and 
SMC gas-to-dust ratios, which we discuss in the 
following section. This plot has been constructed 
in a number of previous works (e.g. Galama & Wi- 
jers 2001; Stratta et. al. 2004; Kami et al. 2006; 
Schady et al. in preparation) but here we show the 
observed distribution of E(B — V) and Nn for the 
first time derived simultaneously from a fit to X- 
ray, optical and nIR data. We find an excess in ab- 
sorption above the Galactic values particularly sig- 
nificant in two sources: GRBs 000926 (E(B — V) 
only) and 010222 (Figure 4), whilst no significant 
intrinsic absorption is necessary in GRBs 970228 
and 990510. The cooling break can be located in 
three of the afterglows: GRBs 990123, 990510 and 
010222 and to all other SEDs a single power law is 
an adequate fit. Details are given below for each 
individual afterglow. 

4.1. Notes on individual sources 

4.1.1. GRB 970228 

No significant absorption is measured for 
GRB 970228. We find only a single power law is 
required and it is possible to pin down the power 
law slope relatively well (y 2 /dof = 11.6/14). 

4.1.2. GRB 970508 

A well-defined temporal decay slope starting 1.9 
days after trigger (Zeh et al. 2006). Preceding this 
time there is an increase in flux followed by an ap- 
parent flattening. The time of the X-ray observa- 
tion log midpoint occurs a little before the 1.9-day 
break at 1.679 days, so we have extrapolated the 
optical data post 1.9 days back to 1.9 days and 
then assume that the evolution of the lightcurve 
is flat back to 1.679 days after trigger. To allow 
for a different behaviour before 1.9 days we in- 
clude a constant value offset in the model between 
the optical and X-ray data, which we both fix at 
1.0 and leave as a free parameter. The improve- 
ment in the fits when the offset is a free parameter 
is somewhat marginal. An F-test indicates that 
the free-parameter model is better at a 98% level, 
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Fig. 3. — Intrinsic absorption in optical/nIR ( E(B — V )) and X-rays (log Nu) measured for the GRB sample. 
We compare the measurements with three different optical extinction laws overlaid with solid curves: Galactic 
(top panel, Predehl & Schmitt 1995), LMC (middle panel, Koornneef 1982; see also Fitzpatrick 1985) and 
SMC (lower panel, Martin et al. 1989). Appropriate metallicities are adopted for LMC (1/3 Zq) and SMC 
(1/8 Zq) calculations (diamonds), and stars mark the centroids of the Solar metallicity fits. For GRB 000926 
the data were too sparse to fit for Nu, so we plot the E(B — V) range at log Nu = 17.0 for clarity. Error 
bars are 90 per cent confidence. 
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Energy (keV) 

Fig. 4. — Spectral energy distributions of GRBs 
000926 (left) and 010222 (right) with best fitting 
models overlaid in red. Overlaid in blue with a 
dashed line is the unabsorbed source flux, demon- 
strating that for these two sources extinction sig- 
nificantly affects the observed optical to X-ray 
emission, in complete contrast to GRBs 970228 
and 990510 where extinction in both optical and 
X-ray regimes are negligible. 


with the offset increasing to 3. Given the uncer- 
tain extrapolation, however, we use the fits with 
the offset as a free parameter in our further anal- 
ysis. The best fitting model for GRB 970508 is a 
single power law with relatively low intrinsic X-ray 
absorption at the level of Nn r 10 21 cm 2 . Kann 
et al. (2006) found a best fit with MW-like dust, 
but as we measure no significant optical extinction 
we cannot distinguish between different extinction 
laws. 

4.1.3. GRB 971214 

GRB 971214 is the highest redshift source in the 
sample at z = 3.418, and we note that the source 
was faint, particularly as seen by LECS, crucial for 
the low X-ray energies. We measure an optical-to- 
X-ray spectral index a ox of 0.6 and find that a 
single power law is an acceptable fit to these data. 
The intrinsic X-ray absorption appears to be ex- 
tremely large whilst the optical extinction is mod- 
erately large, but we note that the high redshift of 
the GRB makes measurement of the X-ray extinc- 
tion more difficult. The curvature in the optical 
part of the spectrum has been previously inter- 
preted as a cooling break (Wijers & Galama 1999) 
and as extinction by either SMC-like extinction 
(Stratta et al. 2004) or by a presently unknown, 
more complex extinction law (Halpern et al. 1998; 
Ramaprakash et al. 1998). Whilst SMC extinction 
is the best fitting law of the 3 used here, it is not 
sufficient to accurately reproduce the shape of the 
optical SED, reflected in the large errors on Nr 
and E(B - V). 

4 .1. 4 . GRB 980329 

In the absence of an accurate redshift determi- 
nation for GRB 980329 we adopt the photometric 
redshift of 2 = 3.6 (Jaunsen et al. 2003), hence 
all results must be taken with caution. The most 
striking feature of this SED is the apparent flux 
deficit in the R band, which is present even after 
correction for Galactic absorption and the proba- 
ble high redshift. We tested for the possibility that 
the R-band flux deficit is due to the 2175A fea- 
ture in the Milky Way extinction curve by fitting 
the pl+MW model whilst leaving the redshift as a 
free parameter lying between 2 = 1.2 and z = 4.2 
(Jaunsen et al. 2003) (we note that at z = 3.6 the 
2175A bump would lie between the I and J bands 
if extinction is MW-like, which is not observed). 
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We obtain a fit with y 2 /dof = 32.7/26 and z = 
1.95lo3i- The other best-fitting parameters, as 
listed in Table 5 for standard 'fits, do not differ 
greatly from those found in the z = 3.6 fit and 
the optical curvature is no better fitted. We also 
tested for the possibility that the R-band deficit 
is caused by a break in the power law, by allow- 
ing the break energy to reside within the optical 
regime (and adopting z = 3.6 with SMC extinc- 
tion). In this fit the break energy could not be 
well constrained, and the power law slopes are not 
well described by the difference of 0.5 as expected 
for a cooling break (the second slope is steeper), 
hence we rule out this possibility. However, it is 
also possible that we are seeing an I band excess 
rather than an R band deficit. Two late time I 
band points taken at 1-10 days since burst (Yost 
et al. 2002) appear not to agree with all of the / 
band data used here (Reichart et al. 1999), which 
could be the result of overestimation of the early I 
magnitudes, underestimation of the later I magni- 
tudes or the occurrence of color evolution. Among 
the models applied to the whole GRB sample here, 
a single power law unabsorbed in the X-rays is the 
best fit, with a moderate E(B — V) of 0.178 and 
a slight preference for SMC extinction. 

4.1.5. GRB 980519 

We caution that all the results are based on a 
redshift estimate equal to the mean of the sam- 
ple spectroscopic redshifts of z = 1.58. A few 
attempts to constrain the redshift have not been 
very precise, amounting to 0.5< z <3.6 or 1.5< 
z <3.6 (lower limits from Jaunsen et al. 2003 and 
upper limit from the fact that we detect a U band 
counterpart). We find that a single power law with 
a small E(B — V) and X-ray absorption consistent 
with zero but with a large error is sufficient to 
model this afterglow. We note that the Galactic 
X-ray extinction towards GRB 980519 is the high- 
est for this sample which, together with the lack 
of known redshift hampers a good measurement of 
A r n for this source. 

4.1.6. GRB980703 

A single power law provides an acceptable fit to 
the spectrum, when absorbed by a large amount at 
both X-ray and optical/UV wavelengths. In this 
afterglow we measure the largest E{B — V) value 
in the sample of 0.29±0.05 or Ay~0.85 at 1.33 


days after burst assuming SMC extinction (which 
is marginally preferred). This is consistent with 
the value found by Bloom et al. (1998’o) of 0.9 ± 
0.2, 5.3 days after the burst. 

A great deal of work has been done on the host 
galaxy properties of this burst because the host is 
bright, with optical extinction measurements by 5 
groups. There is a hint that the optical extinction 
may be decreasing with time (see e.g. Holland 
et al. 2001) since measurements of Ay at differ- 
ent times are inconsistent: 2.2 at 0.9 days, 

Castro-Tirado et al. 1999; 1.5 at 1.2 days, 

Vreeswijk et al. 1999; Ay~0.3 at 4.4 days, Djor- 
govski et. al. 1998; -~0.9 at at 5.3 days, Bloom et 
al. 1998b, and there appears to be a discrepancy 
between the measured optical spectral and tem- 
poral slopes when assuming Ay is constant. How- 
ever, we note that the optical spectral slope was 
taken to be (3qa = -2.71i0.12 from Vreeswijk et al. 
(1999), and in this study we obtain a lower value 
of Poa = -0.92i0.03 which would be completely 
consistent with /?ox = (li2a)/3 = 0.9 using a = 
0.85 from Zeh et al. (2006, noting their a value 
has a very large associated error - Table 1). We 
use the Vreeswijk et al. (1999) optical data here 
and scale it from 1.2 days to 1.33 days after trig- 
ger. Combining the optical and X-ray data when 
fitting provides us with a different estimate for the 
extinction than was obtained by Vreeswijk et al. 
for the optical data alone. We note that we have 
made only a minimal extrapolation of the original 
optical data used in this analysis, from ~1.2 days 
to 1.3 days after trigger. 

4.1.7. GRB 990123 

GRB 990123 does not have significant X-ray ab- 
sorption above the Galactic value, and the optical 
extinction is consistent with zero. We set an upper 
limit to the latter which is ten times lower than 
the value found by Savaglio et al. (2003) from fits 
to Zn II and H I in the optical spectrum. A single 
power law fit to both optical and X-rays results in 
a spectral slope of /3=0.61±0.01 at 1.24 days since 
burst, comparable to the /?ox=0.67±0.02 at about 
the same time since the burst found by Galama et 
al. (1999). The latter authors note that the cool- 
ing break must lie at or above X-ray frequencies 
at that time. However, we find an improved fit 
with a broken power law, constraining the cool- 
ing break to 0.15 < v c < 2.4 keV, within the X- 
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ray spectrum, also found by Stratta et al. (2004), 
Corsi et al. (2005) and Maiorano et. al. (2005). 
We note, however, that the X-ray spectrum com- 
prises flux accumulated between 0.2 and 2.6 days 
since burst, including the proposed jet break time 
of 2 days (Table 2). We tested for the possibil- 
ity that the offset between optical and X-ray data 
was incorrect, but this made very little difference 
to our overall goodness of fit. It is possible that the 
cooling break has entered the X-ray band during 
these observations, since this break is expected to 
decrease in frequency with time. If this is the cool- 
ing break, the spectral slopes above and below the 
break differ by the expected factor of 0.5 when left 
free. We therefore took only the data from the first 
X-ray observation, before the jet break (Zeh et 
al. 2006), together with scaled optical/nIR data, 
which has a log midpoint of ~0.7 days since trig- 
ger. We repeated all our fits. We find acceptable 
fits again only with the broken power law model, 
with the break energy lying somewhere between 
the optical and ~2 keV (preferring central values 
just below X-ray frequencies, in line with our X- 
ray only analysis in which a single power law is 
a good fit to all spectra). We note that an ex- 
cess of flux at high energies (as seen by the Bep- 
poSAX PDS instrument not used here) is reported 
by Corsi et al. (2005). They attribute this to 
an Inverse Compton component (though we note 
that this remains inconsistent with the radio data, 
Kulkarni et al. 1999b). Since the X-ray spectrum 
is adequately fit with a single power law we as- 
sume that the tail of any such component is not 
significant below 10 keV. 

11.8. GRB 990510 

Fits with Galactic, LMC- and SMC-like extinc- 
tions show that E(B — V) is very low in this source 
and we can only determine upper limits. The low 
amount of extinction makes it impossible to pin 
down the extinction curve shape, hence a fit us- 
ing Galactic-like extinction is sufficient here and 
LMC- and SMC-like extinctions give similar re- 
sults. There is considerable improvement in the 
X 2 when allowing for a break in the power law, 
noted by previous authors, which we find is located 
at 0.016-0.020 keV at ~1.06 days since burst (of 
the order of the value of ~0.029 keV at ~1 day 
since burst given by Pian et al. 2001). The slope 
change is as expected for a cooling break in the 


slow cooling regime when leaving both power law 
slopes free. Our measured spectral slope in the 
optical regime agrees well with that measured'in 
the Very Large Telescope (VLT) spectra of ft = 
0.6±0.1 by Vreeswijk et al. (2001). 

We note that Kuulkers et al. (2000) analysed 
these X-ray data in several time bins and found 
no spectral evolution, hence the cooling break re- 
mains outside the X-ray frequencies during the ob- 
servations. We constrain the optical extinction to 
be E(B — V)< 0.01, and the X-ray equivalent hy- 
drogen column to be < 0.15xl0 22 cm -2 for SMC 
extinction and < 0.87xl0 22 cm -2 for MW extinc- 
tion. From the X-ray spectrum alone a higher col- 
umn of Nh = 2.1 ± O.GxlO 21 cm -2 was measured 
by Kuulkers et al. (2000). Optical spectra have 
provided a lower limit on the amount of neutral 
hydrogen towards GRB 990510 of log JV(H I) > 
19.7 cm -2 (Vreeswijk et al. 2001). These authors 
obtain an approximate estimate for the metallic- 
ity from the optical spectra using Fe/H, and find 
12+ log [Fe/H] = -1.5 ± 0.5 or 0.01-0.1 times the 
Solar value. This range approximately covers the 
metallicity of both the LMC (0.33, Pei 1992) and 
the SMC (0.125, Pei 1992) so we also fitted the 
data with these two metallicities. As there is no 
substantial absorption observed, the fits do not 
change significantly. 

4-1.9. GRB 000926 

The optical and IR SED of this burst is very 
well sampled, but the X-ray afterglow was very 
faint at the time of observation. For this rea- 
son these X-ray data have been grouped to have a 
maximum of only 10 counts per bin, and strictly 
speaking this means Gaussian statistics should be 
treated with caution. However, we do use the x 2 
statistic as a goodness of fit for comparison with 
the rest of the sample. Given the poor quality of 
the X-ray data, and the lack of sufficient counts 
in low X-ray energy bins, we have fixed the X- 
ray column density at zero and fit only for the 
optical extinction. We find a large amount of in- 
trinsic extinction is necessary to describe the flux 
deficit with respect to a single power law, consis- 
tent with the value found by most other studies 
(Stratta et al. 2004; Fynbo et al. 2001; Harrison 
et al. 2001; Price et al. 2001). A large extinc- 
tion was also derived from an optical spectrum 
by Savaglio et al. (2003), and the Ay found in 
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that study is approximately twice the value found 
here. However, Savaglio et al. use the spectral 
line measurements to first fit for the depletion pat- 
tern and then infer an extinction. Harrison et al. 
(2001) interpreted the optical flux deficit as in- 
dicating that a significant fraction of the X-ray 
flux was in fact Inverse Compton emission - later 
also suggested for GRB 990123 (see above). This 
is the only afterglow for which an LMC extinc- 
tion law is (marginally) preferred. Fynbo et al. 
(2001) report a tentative H I column density mea- 
surement of A(H I) ~2xl0 21 cm -2 which leads to 
a relatively high metallicity with [Zn/H] = -0.13. 
This metallicity is between the LMC and the MW 
values, which may explain the preference here for 
LMC-like extinction if no 2175A bump is present. 

4.1.10. GRB 010222 

A good dataset for GRB 010222 allows the spec- 
tral properties to be well constrained. We find 
that whilst a single power law is a reasonable fit 
to these data, a broken power law significantly im- 
proves the fit. The break energy lies around 0.01 
keV, above the frequency of the last optical band 
in our SED. Optical extinction is clearly non-zero, 
with E(B — V) = 0.06±0.02, consistent with that 
found by Lee et al. (2001) for an SMC extinction 
law but about three times lower than that inferred 
from the spectral lines by Savaglio et al. (2003). 
X-ray absorption is also required with an effective 
hydrogen column of A h = 1.15+®'^ x 10 22 cm -2 . 

Panaitescu & Kumar (2002) in their fits to 
the multiwavelength afterglow of this source find 
significant reddening of the optical spectrum of 
Ay — 0.21 with an SMC extinction curve, which 
is consistent with our value. But they find a 
large fitting error and attribute this to 8 outlying 
points suggesting that either some reported obser- 
vations have underestimated uncertainties or there 
are short timescale fluctuations in the afterglow 
emission (Cowsik et al. 2001). Their jet model re- 
quires the cooling break to pass through the X-ray 
band at about 1 day, which they find incompat- 
ible with their observations. Our analysis places 
the cooling break at optical/UV wavelengths at 
1.51 days since burst. 


4.2. Comparison with previous studies 

We can compare our results directly with those 
of previous studies of samples overlapping with 
this BeppoSAX subsample. In general, we are find- 
ing similar central values for extinction as all pre- 
vious studies, and are improving upon the uncer- 
tainties, fitting all afterglows in the same consis- 
tent manner allowing for direct comparison. 

Galama & Wijers (2001) performed the first 
systematic study of line-of-sight Ah and E(B — V) 
with a sample of 8 afterglows, consisting of all but 
the two most recent bursts in our sample. From 
fits to the X-ray spectra they found intrinsic Ah 
amounting to 10 22 -10 23 cm -2 , ruling out the pos- 
sibility that some hosts have no X-ray column at 
all. They noted that these values lie in the range of 
Galactic giant molecular clouds (estimating cloud 
sizes of 10-30 pc) - a conclusion also recently ar- 
rived at when including Swift bursts (Campana et 
al. 2006a) and when measuring H I from damped 
Lyman alpha absorption in GRB optical afterglow 
spectra (Jakobsson et al. 2006a). They used a 
simple extinction law Ay oc v with a smoothly 
broken power law. Comparing their dust to gas 
ratios with that of the Milky Way they obtain an 
optical extinction 10-100 times smaller than ex- 
pected. Their finding of generally low Ay was at- 
tributed to previously predicted dust destruction 
by the GRB. Of the bursts common to both their 
and our samples and having a known redshift we 
obtain consistent extinction values within the 20- 
uncertainties in all cases except for GRB 971214 
for which the two optical extinctions are only con- 
sistent at the 3o- level, our central value lying four 
times lower than the Galama & Wijers measure- 
ment. We obtain smaller uncertainties in our ex- 
tinction measurements in all cases except for the 
X-ray column of GRB 980703. 

Stratta et al. (2004) have also measured Ah 
and Ay r in the optical and X-ray data seperately, 
and later plotted the combined data in flux space 
(after assuming the X-ray model to be correct) in 
order to judge the position of the cooling break. 
To derive optical spectral slopes, Stratta et al. 
first adopt the p-value derived from the X-ray fits 
(the input electron energy index, see Paper II for 
our fits for this parameter), use this to fix the op- 
tical spectral index and fit for Ay . Our results for 
E(B — V) with SMG-like extinction are consistent 
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in all cases to within the 90 % confidence limits, 
adopting the cooling break positions found from 
this work if known (the values for GRB 970508 
are only consistent if v a > u c ). Even for the two 
GRBs for which we have assumed different red- 
shifts we find consistency in extinction estimates 
(Stratta et al. adopt z = 1 for both 980329 and 
980519). In 7 of the 10 cases we derive a better 
constrained value or upper limit to the extinction. 
The X-ray absorbing columns we measure are also 
generally consistent with those found by Stratta 
et al. in a fit to the X-ray data only (these are 
of course the same data used in our study, except 
for differences in number of observations combined 
for some sources, and we include GRB 000926 as 
well). Our method obtains improved estimates for 
Nh for 7 of the 9 sources common to both stud- 
ies. For GRB 990510 Stratta et al. find a value 
higher than our derived 1V H using Solar metallic- 
ity, and for GRB 980329 our result for 1V H is less 
accurate at the 90 % confidence level, though we 
note the different redshifts assumed hence direct 
comparison is not possible. 

Kann et al. (2006) fitted only the optical SEDs 
for a sample of pre-Swift bursts including 8 anal- 
ysed here (not including GRBs 970228 and 980329 
and using z — 1.5 for GRB 980519). We find 
similar central values of Ay and improve upon 
their mean error by a reduction of 5-10 %. Val- 
ues disagree at the 90 % confidence level only for 
GRB 970508, where we find an upper limit to the 
extinction which is 2.4 times lower than the lower 
limit of Kann et al. but would be consistent with 
their estimate at the 3cr level. 

5. Discussion 

For half the afterglows the best-fitting model to 
the SED includes SMC-like extinction (as opposed 
to LMC or MW) and in one case LMC-like extinc- 
tion. In no cases is there a preference for MW-like 
extinction. We are sensitive to the 2175A bump 
(MW) in the redshift range z = 0.46-9.9, cover- 
ing all our selected GRBs, but clearly we do not 
detect any such feature. We find a wide spread in 
central values for the gas-to-dust. ratios, and for 4 
GRBs the gas-to-dust ratios are formally inconsis- 
tent with MW, LMC and SMC values at the 90 
% confidence limit assuming the SMC metallicity 
(Figure 5). In these 4 cases the ratio is several 


orders of magnitude higher than the SMC value 
of 4.4±0.7xl0 22 cm -2 mag" 1 (Koornneef 1982; 
Bouchet et al. 1985) and must mean 'that ei- 
ther gas-to-dust ratios in galaxies can span a far 
larger range than thought from the study of lo- 
cal galaxies, or the ratios are disproportionate in 
GRB hosts because the dust is destroyed by some 
mechanisms (likely the GRB jet), or that the lines 
of sight we probe through GRBs tend to be very 
gas-rich or dust.-poor compared with random lines 
of sight through galaxies. Finally, a dust grain 
size distribution which is markedly different than 
considered here may also affect these ratios. 

In fact, a recent study has shown that for the 
LMC the former is true. A recent observation of 
four core-collapse supernova remnants (SNRs) in 
the LMC with Spitzer has shown IR emission as- 
sociated with the supernova blastwave (Williams 
et al. 2006). This is interpreted as dust with an 
LMC-like grain size distribution which has been 
collisionally heated by the X-ray emitting plasma. 
The observations require that some fraction of the 
small dust grains has been destroyed by sputtering 
by high energy ions in fast shocks. Dust destruc- 
tion is known to occur in SNR shock fronts (Jones 
2004), and we will return to this issue later in the 
section. The derived gas-to-dust ratios are sev- 
eral times higher than the LMC ratio, as we see in 
the line-of-sight measurements of GRBs and has 
been observed in other types of supernova (e.g. 
Borkowski et al. 2006), the cause of which is not 
known. 

We measure a large amount of intrinsic absorp- 
tion in some of the sample (Figure 4), and can 
state that absorption is insignificant in others, as 
seen for example in the contrast between GRBs 
010222 and 990510. We have tested for the possi- 
bility that the Galactic column density in the GRB 
direction is affecting our intrinsic column deriva- 
tions b 3 r plotting JV H ,int vs. JVh, G al, and found 
no correlation. We assumed therefore that all the 
measured absorption lies within the host galaxy. 
We note that a few of these afterglows have spec- 
tra: for GRB 990510 no intervening systems were 
clearly identified, an intervening system is mea- 
sured at 168 km s -1 from the host redshift for 
000926 (assumed to lie within the host galaxy, 
Castro et al. 2003), 2 intervening systems are 
found for GRB 010222 (Jha et al. 2001; Mirabal 
et al. 2002) and at the host redshift there are two 
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components separated by 119 km s _1 (Mirabal et 
al. 2002). None of these intervening systems are 
close enough to us as the observe!' to significantly 
affect measurement of the intrinsic host extinction. 

5.1. Approaches to measuring absorption 
in the host galaxies 

In this study we provided a thorough, uniform 
study of both optical and X-ray extinction along 
the lines of sight towards a sample of 10 GRB af- 
terglows. The well known spectral shape and rel- 
ative brightness of the afterglow emission make 
GRB afterglows a powerful line-of-sight probe of 
high redshift extinction. This is one of several ap- 
proaches to measuring absorption in GRB host 
galaxies. One can globally divide the studies of 
extinction in the field of GRBs in two categories: 
line-of-sight extinction studies and studies of ex- 
tinction of the integrated host galaxy or parts 
thereof. 

5.1.1. Line-of-sight studies 

Line-of-sight studies generally involve fitting 
the afterglow spectral energy distributions in opti- 
cal and/or X-rays with template extinction mod- 
els (i.e. MW, LMC, SMC or more parametrized 
models) as we have done here. The standard 
conversions between X-ray extinction and opti- 
cal extinction for the Milky Way and the two 
Magellanic Clouds are generally in disagreement 
with the column densities measured through this 
method, but these skewed gas-to-dust ratios are 
also being found in other astrophysical situations 
as discussed above for SNRs in the LMC, and the 
destruction of dust can go some way to alleviating 
the mismatch. 

A further way to probe line-of-sight extinction 
properties is through optical spectroscopy of the 
afterglow. In this case the careful measurement 
of column densities of heavy elements can be used 
to study the dust depletion pattern along the line 
of sight (e.g. Savaglio & Fall 2004). The mea- 
sured metal column densities in combination with 
the best-fitting depletion pattern and the empiri- 
cally determined conversion between Ay and the 
dust column, can provide a prediction of the dust 
extinction along the line of sight to a GRB (for 
a detailed explanation see e.g. Savaglio & Fall 
2004). Savaglio & Fall (2004) show that the ex- 


tinction derived from the dust depletion method 
is significantly higher than the value derived from 
direct fitting to the continuum of the afterglow 
spectrum, e.g. for GRB 020813 they find an over- 
estimation by at least five times. We show that 
this conclusion holds when fitting the afterglow 
continuum emission over a much larger wavelength 
range and can quantify that overestimation factor 
for GRBs 990123, 000926 and 010222 as approx- 
imately > 11, > 2 and > 3 times overestimated 
respectively. The reasons for this apparent dis- 
crepancy may be two-fold: firstly the fitted ex- 
tinction profiles to the afterglow SED are likely 
poor approximations to the true extinction pro- 
file, and secondly the host dust depletion chem- 
istry may well differ from the Milky Way chem- 
istry. In addition, the GRB or afterglow may pref- 
erentially destroy small dust grains, skewing the 
extinction profile towards larger grains, resulting 
in a ’’grey” extinction curve. This would alter the 
derived extinction from SED fitting, and possibly 
bring estimates from dust depletion methods and 
SED fitting closer together. Whist a grey extinc- 
tion curve was the best fitting extinction curve to 
GRB 020405 (Stratta et al. 2005), we note that 
grey extinction curves have been fitted to samples 
of afterglow SEDs with no conclusive improvement 
in fit (e.g. Stratta et al. 2004). 

5.1.2. Integrated host galaxy studies 

One can also study the extinction properties of 
host galaxy as a whole, and there are again sev- 
eral methods to do this. Whilst it has been shown 
that GRBs occur in starforming regions in the host 
(Bloom et al. 2002; Fruchter et al. 2006), many 
host galaxies are small and mixing timescales may 
be short, enabling global properties to be mea- 
sured. The host and afterglow have similar, mod- 
erate reddening in GRB 000418 (Gorosabel et al. 
2003a) which is taken as evidence that the ISM is 
well mixed. But more extreme values of reddening 
are also seen, such as the extremely red afterglow 
and host of GRB 0301 15, in which the host is an 
Extremely Red Object (ERO, Levan et al. 2006). 

One of the most common methods of integrated 
host galaxy studies is fitting of the broadband op- 
tical and near-infrared SEDs of the hosts them- 
selves (e.g. Christensen et al. 2004). Galaxy tem- 
plates can be fit to the data, using photometric 
redshift programmes such as HyperZ (Bolzonella 
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et al. 2000) providing values for the photometric 
redshift, the age of the dominant stellar popula- 
tion and the extinction, by fitting a series of galaxy 
templates. The extinction measured this way is 
the extinction by the ISM of the galaxy on the stel- 
lar light,, E(B — V) s , in which the geometry of the 
dust in the galaxy can play an important role (e.g. 
dustlanes as opposed to an homogeneous dust dis- 
tribution). A study of a large sample of hosts 
has been performed through HyperZ template fit- 
ting by Christensen et al. (2004), who find that 
GRB hosts generally exhibit little extinction, and 
have young stellar populations. The dependence 
on metallicity and assumed initial mass function 
is small (Gorosabel et al. 2003a, b; Christensen et 
al. 2004). One of the difficulties faced here is that 
the galaxies are often very faint. 

5.1.3. Emission-line spectroscopy 

In low redshift (z < 1) cases, an optical spec- 
trum of the host galaxy can be taken, typically 
showing several nebular and Balmer emission 
lines. The Balmer lines can be used to derive 
values for the reddening by calculating their devi- 
ation from case B recombination values expected 
in typical starforming region conditions (Oster- 
brock 1989). The derived reddening E(B — V) g is 
the reddening of the ionized gas in the source, i.e. 
the dominant starforming region(s) producing the 
Balmer emission lines. In general the reddening 
found from the Balmer decrement is low (see e.g. 
Prochaska et al. 2004), to very low (Wiersema et 
al. 2006). 

The two host galaxy extinction estimates 
E(B — V) g and E(B — V) s may be correlated 
for most galaxies (see Calzetti 2001 for a review). 
The increasing data volume on nearby GRB host 
galaxies will allow a test of these correlations, 
providing further insight into effective GRB host 
galaxy extinction curves. 

On rare occasions it is possible to obtain high- 
resolution spectra of an afterglow that also shows 
host galaxy emission lines, allowing one to obtain 
a simultaneous view of the extinction along the 
line of sight and of the Balmer decrement. In 
the case of GRB 060218, both absorption lines and 
emission lines are detected at high resolution using 
the UVES spectrograph on the VLT (Wiersema et 
al. 2006). The spectrum shows asymmetric emis- 
sion lines that are well fitted with two Gaussians 


separated by 22 km s _1 . The same two veloc- 
ity components can be seen in absorption in Ca 
II and Na I, and have different, chemical proper- 
ties. These two systems can be interpreted as two 
seperate starforming regions, through which the 
light of the afterglow shines. A broadband mea- 
sure of the extinction either from the afterglow or 
from template fitting of the host would not have 
been able to seperate out the contributions of the 
two individual systems. 

5.1.4- The longer wavelengths 

Yet another way to detect the presence of dust, 
is the detection of GRB host galaxies in the far- 
infrared or sub-nun (see e.g. Barnard et al. 2003; 
Tanvir et al. 2004), where the UV radiation from 
massive stars is reprocessed by dust and reradi- 
ated in the far infrared. Detection of hosts in the 
far infrared as well as optical can severely con- 
strain their SED (Le Floc’h et al. 2006), providing 
estimates for the unobscured star formation rate 
of GRB hosts. In a few cases values up to hun- 
dreds of solar masses per year have been reported, 
while optical indicators give much lower values, 
indicating a lot of dust,- obscured star formation. 
A different probe of unobscured star formation is 
the radio continuum flux, which is thought to be 
formed by synchrotron emission from accelerated 
electrons in supernova remnants and by free-free 
emission from H II regions (Condon 1992). Berger 
et al. (2003) performed a survey of host galaxies at 
radio and sub-mm wavelengths, finding that a sig- 
nificant fraction of GRB host galaxies have a much 
higher radio-derived star formation rate than opti- 
cal methods indicate, pointing again to significant 
dust ext.incted star formation but based on only a 
few detections. 

5.1.5. Our method compared 

All methods considered, SED fitting of the af- 
terglow as described and carried out in this pa- 
per is the most broadband view of line-of-sight, 
properties. Given the very well known underly- 
ing continuum shape, the extinction curve can in 
principle be extremely well modelled using this 
method. As discussed in Section 1, the effects of 
dust extinction in the optical are better measured 
at higher redshifts whereas the X-ray absorption 
is best measured at lower redshifts, so this method 
is most reliable for some middle range of redshifts 
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centred around z = 2 which approximately corre- 
sponds to the current Swift median redshift (e.g. 
Jakobsson et al. 2006b). As with any line-of- 
sight method, the measured columns may not be 
representative of the host galaxy as a whole, so 
comparison with the integrated host galaxy meth- 
ods is important. The general drawbacks of this 
method lie in the inability to disentangle the loca- 
tions of absorbers, both intervening (between us 
and the host) and within the host itself. Inter- 
vening absorbers may be identified using afterglow 
spectroscopy (Ellison et al. 2006). Distribution of 
absorbers within the host may be tackled from the 
point of view of searches for variability in the ab- 
sorption which may indicate destruction of dust 
or ionisation by the GRB and/or afterglow, and 
comparison of the measured columns with known 
absorbers such as Molecular clouds. Campana et 
al. (2006a) find that the X-ray absorptions of 17 
bursts (including most of this sample) are consis- 
tent with the GRBs lying along random sightlines 
towards Galactic-like Molecular clouds. Our sam- 
ple is too small to make such a statistical compar- 
ison, but we note that the Campana et al. sample 
includes the results of the Stratta et al. (2004) 
study of the GRBs presented here for which we 
obtain largely consistent extinction values. 

Dust destruction is of course impossible to mea- 
sure in a single SED, but may be seen with 
multiple epochs of data. Destruction of grains 
by for example the afterglow UV/X-ray emission 
would in principle be observable, because deple- 
tion indicative line ratios would vary on observable 
timescales as metals are released from grains into 
the ISM. Variable line emission has been searched 
for in afterglows with multiple optical spectra and 
has been seen in only two cases, GRB 020813, 
in which the Fe II A2396 transition equivalent 
widths decreases by at least a factor of five over 
16 hours (Dessauges-Zavadsky et al. 2006) and 
GRB 060418 in which several transitions of Fe II 
and Ni II are seen to vary in observations cover- 
ing 11 to 71 minutes post-burst (Vreeswijk et al. 
2006). This line emission is thought to arise at 50- 
100 pc from the GRB site, possibly within range 
of the GRB ionising flux. A variable Ay has only 
been suggested for GRB 980703 (described in Sec- 
tion 4.1.6). There is evidence for time variable X- 
ray absorption in a small number of GRBs ranging 
from tentative to moderately strong, in particular 


for GRBs 011121 (Piroet al. 2005), 050730 (Star- 
ling et al. 2005) and 050904 (Boer et al. 2006; 
Campana et al. 2006b; Gendre et al. 2006) and 
060729 (Grupe et al. 2006). This implies ionisa- 
tion of the line-of-sight gas, probably by the high 
energy GRB jet (though we note that column den- 
sity changes and early-time spectral evolution can 
have similar effects on the spectrum). 

Dust, destruction can, however, only affect the 
immediate environment around a burst, (perhaps 
out to 10 pc, Waxman & Draine 2000), and not all 
the dust in that region need be destroyed since this 
depends on the efficiency of the mechanism and 
on the dust grain size distribution. Kann et al. 
(2006) noted a possible tentative correlation be- 
tween Ay and star formation rate as measured via 
host galaxy emission in the sub-mm (for a sample 
of 7 GRBs) . If this turns out to hold when tested 
against larger samples it would demonstrate that a 
substantial amount of the dust we see comes from 
star forming regions located throughout the whole 
host galaxy rather than very close to the burst 
where it could be destroyed. However, this corre- 
lation is unlikely to be real given the very low res- 
olution of the sub-mm observations which make it 
difficult to differentiate between host galaxy emis- 
sion and field galaxy emission. This differentiation 
can be made with Spitzer and lower star forma- 
tion rates in comparison with the sub-mm derived 
values have been found for some hosts including 
those of GRBs 980703 and 010222 (Le Floc’h et 
al. 2006). 

6. Conclusions 

Here we have demonstrated the advantages of 
simultaneous fitting of broadband data using a 
subsample of BeppoSAX GRB afterglows. In no 
cases is a MW-like extinction preferred, when test- 
ing MW, LMC and SMC extinction laws. The 
2175A bump would in principle be detectable in 
all these afterglows, but is not present in these 
data. An SMC-like gas-to-dust ratio (or lower 
value) can be ruled out for 4 of the hosts anal- 
ysed here (assuming SMC metallicity and extinc- 
tion law) whilst the remainder of the sample have 
too large an error to discriminate. 

We discuss the various methods employed to 
derive host galaxy extinctions, and compare our 
results with previous works. We find that this 
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Fig. 5. — Gas-to-dust ratios, Nu/E(B — V), de- 
rived from best fits to the SEDs assuming SMC 
metallicity, plotted for each sample GRB in date 
order from left to right (GRB 970228 first), exclud- 
ing GRB 000926 where we did not fit for X-ray 
column density. We include the ratios for 3 GRBs 
taken from the literature for comparison (square 
data points): GRBs 000301c (Jensen et al. 2001), 
000926 (Fynbo et al. 2001) and the lower limit 
on 020124 (Hjorth et al. 2003). The solid, dot- 
ted and dot-dashed lines show the measured gas- 
to-dust ratios and their errors for the Milky Way 
(Diplas & Savage 1994), LMC and SMC (Koorn- 
neef 1982; Bouchet et al. 1985). Error bars are at 
the 90 % confidence level. 


method provides similar central values of E(B—V) 
and Nh to previous works in which extinction or 
absorption is determined through afterglow con- 
tinuum fitting, and in the majority of cases we 
obtain tighter constraints. We confirm that, with 
respect to continuum fitting methods such as this, 
optical extinction is overestimated with the deple- 
tion pattern method, and quantify this for a small 
number of cases. 

Swift, robotic telescopes and Rapid Response 
Mode on large telescopes such as the William Her- 
schel Telescope on La Palma and the Very Large 
Telescopes in Chile now allow earlier and higher 
quality data to be obtained, which will help im- 
mensely in discriminating between the different 
extinction laws at work in the host galaxies. 
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